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Substructures formation induced by gravitational tides? 



F. Renaudi'2, Ch. Theis^, T. Naab^ and C. M. Boily^ 

Abstract. Physics lectures always refer to the tides as a disruptive effect. 
However, tides can also be compressive. When the potential of two galaxies 
overlap, as happens during a merger, fully compressive tides can develop and 
have a strong impact on the dynamics of substructures such as star clusters or 
tidal dwarf galaxies. Using A^-body simulations of a large set of mergers, we 
noticed the importance of these tidal modes at cluster scale. With a model of 
the Antennae galaxies, we conclude that the positions and timescales of these 
tidal modes match the actual distribution of young clusters. A detailed study 
of the statistics of the compressive tides shows a stunning correlation between 
this purely gravitational effect and the observed properties of the star clusters. 
In this contribution, we introduce the concept of compressive tide and show its 
relevance in the special case of the Antennae galaxies. We extend our conclusions 
to a broad range of parameters and discuss their implications on several critical 
points such as the infant mortality, multiple star formation epochs in clusters or 
the age distribution. 



1 Introduction 



The role of the tides in the evolution of galaxies has been known to be of prime 
im portance since the earliest observations and simulations of such systems (see 
e.g. Sandage 1961; Toomre & Toomr ?1972'; 'To omrdl 19771 : iHibbard fc van GorkomI 
Il99 6: Lai ne et al. 20Q3,; Elmegreen et al. 20071 )7" With increasing resolution in 
both methods, detailed investigations on the structures formed in galactic pairs 
became possible. Among them, the study of star clusters and tidal dwarf 
galaxies (TDGs) gathered a lot of efforts with the closest major merge r as a 



reference: the Antennae ga l axies (see e.^ pj. Barnes _ 198^ : Mirabel et al 



Whitmore Schweizer 1995: 'Meureilll995: Hibbard et al.'' 200ll : lFall et al. 



1992; 



2005 



Mengel et al.ll2005 : .Whitmore et al.ll2007 : Karl et al. 2008l '). Recent HST obser- 



vations of this pair revealed an amazing firework of star formation in the central 
regions. 

Though established, the link between the parsec-scaled structures like the 
young star clusters, and the phenomena covering over tens of kiloparsec, as tides 
do, still misses a quantitative description. \ ye pr opose to partially bridge the 
gap between these two scales. iRenaud et al.l ( 20081 ) noted that fully compressive 
tides develop in interacting galaxies and have characteristic timescales of some ^ 
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10 Myr, matching well those of the formation and early life of star clusters. They 
proposed that this additional gravitational effect would increase the binding 
energy of a young cluster and thus provide a favorable environment for the 
formation of stars. 

Considering these tidal modes as the boundary conditions of a cluster, we 
propose that they help to form and to bind a cluster, and this way, preve nt it 
from early dissolution ( Whitmore et al.ll2007l : Ide Griis Sz Parmentierll2007l ). In 



this contribution, we extend the conclusions of Renaud et al. ( 20081 ) thanks to 



a parameter survey. After a brief description of the method used, we present 
the results of the Antennae galaxies and then, the influence of the pericenter 
distance between the progenitors on the statistics of compressive tidal mode. 
We close with a discussion on the role of the compressive tides on the evolution 
of star clusters. 



2 Numerical method 

The properties of the tides rely on the shape of the gravitational potential. To 
derive them, we use collisionless iV-body simulations of typical galaxy merg- 
er^. The resolution used for this study counts 1,400,000 particles shared among 
the disk-bulge-halo components of the models. The smoothing length and the 
number of particles match the typical cluster scale (40 pc a nd 10^ M(7)). 

The equations of motion are integrated by gyrf alcON ('Dehnen"2002') . The 



tidal tensor (i.e. the second spatial derivative of the potential, see Ren aud et al, 



20081 for details) is computed at the position of a stellar particle and then diago- 
nalized. The eigenvalues determine the strength of the tides along the associated 
eigenvectors. The nature of the tidal field only depends on the sign of the largest 
eigenvalue: when positive, the tides tend to stretch the matter along at least 
one direction; when negative, all the eigenvalues denote a compression and the 
tides are fully compressive. 

By applying this method via efficient algorithms to all the stellar particles, 
we derived statistics of the compressive tidal mode, and its evol ution in time. 



Deta ils of these methods will be presented in a forthcoming paper ([Renaud et al.l 
20091 1. 



3 The Antennae 



Our st arting point is the model of the Antennae galaxies presented in I Renaud et al 
(|2008l ). T he two self-gravita ting equal-mass progenitors have been s et up with 
magalie ( Boilv et al. 200ll ) which gathers an exponential disk, a Hernquist 



([1990) bulge and an isothermal dark matter halo, all scaled to get two SO galax- 
ies. The progenitor have been set on a bound orbit with an eccentricity e ~ 0.96 
and an initial separation of 75 kpc. 

In Fig. [H we show the evolution of the stellar mass fraction in compressive 
mode. Before the first passage {t < 0), the galaxies are isolated and have an 



^ The gas, which represents a small fraction of the mass, is here supposed to follow the stellar 
component and thus will not affect strongly the gravitational potential. 
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Figure 1. Evolution of the mass fraction in compressive mode for our Anten- 
nae model. The equivalent morphology is shown for three snapshots (50 Myr, 
150 Myr and 300 Myr) with column density (in blue) and the compressive 
particles (in red). The two pericenter passages (t = and t = 300 Myr) are 
associated with the two main peaks, rising the mass fraction in compressive 
mode to « 14%. 



intrinsic compressive fraction of ~ 3%. The main tidal events like the creation 
of bridges and tails {0 < t < 100 Myr) and the merger phase (250 Myr < t < 
450 Myr), are clearly visible on this plot, as they are linked to a raise of the mass 
fraction by a factor ~ 5. Note that the compressive zones spread in the nuclei 
area but also in the tails , and in particular n ear the TDGs region (jSchweizer 
197i : iMirabel et aDll992l : iHibbard et alll2005l ). 

Typically, one of our stellar particle enters and leaves compressive regions 
and thus undergoes a series of on/off episodes along its orbit. The age distribu- 
tion of these events follows an exponential law with a timescale of ~ 10 Myr. By 
considering a mean distribution between continuously compressive events and 
on/off episodes, it is possible to retrieve a power-law dN/dr oc t~^'^ age distri- 
bution of the compressive modes. With the rough assumption that star clusters 
would dissolve when the tidal field switches from compressive to extensive, one 



(|2005l . 



can directly compare this distribution with the observations of lFall et al 
dN/dr oc t^^) for the age of the Antennae clusters. 

Because of a growing number of clues suggesting a link between compressive 
tides and young clusters, we decided to extend our conclusions to a larger set 
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Figure 2. Evolution of the mass fraction in compressive mode for different 
pericenter distances d, normalized to the exponential disk scalelength Rd of 
the progenitors. The fraction increases at the pericenter passage(s) when the 
interaction is penetrating. In the other cases, it remains almost fiat, at the 
level of the isolated progenitors (~ 3%). 



of parameters. Among them, we focus here on the influence of the pericenter 
distance. 



4 Influence of the pericenter distance 

In this section, we investigate the role of the distance between the two progeni- 
tors at the first passage. For simplicity, all the other parameters (scalelengths, 
inclination, spin, eccentricity ...) remain as in the Antennae model. We set the 
orbits of the galaxies to cover the full range of distances, from mergers (short 
distance at first passage) to fly-bys. The Fig. [5] plots the evolution of the mass 
fraction in compressive mode (same as Fig. [1]) for different pericenter distances 
d, normalized to the scalelength of the exponential disk i?^ of the progenitors. 

Again, the main steps of the history of the interaction are clearly visi- 
ble. However, we note that the mass fraction in compressive mode decreases 
as the pericenter distance grows. As for the Antennae galaxies {d/Rd = 1.3, 
not shown here), when one progenitor goes deeply into the potential well of the 
other {d/Rd = 2,4), the fraction in compressive mode increases by a factor four 
to five. On the contrary, wider interactions show a milder increase {d/Rd = 8) 
or remain at the level of isolation {d/Rd = 14). In this fashion, if the compres- 
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sive tidal modes are to be linked to an enhanced star formation rate (SFR), we 
would expect a higher rate for close passages, thus mergers. 

^Other studies showed an enhanced SFR, even for fly-bys (see e.g. ldi Matteo et al. 

20071 ). This difference may be interpreted as gas stripped off one galaxy to the 
other, a phenomenon that occurs for all types of interactions. Yet, it is difficult 
to predict an accurate fate of the giant molecular clouds or the young star clus- 
ters, as their response to a compressive non-isotropic non-constant tidal field 
has not been investigated so far. 



5 Compressive tides and multiple populations 



The statistical study of the tidal field revealed global properties that ha ve been 
extended to the general case of galaxy mergers (iRenaud et al.ll2008l . l2009l ^ . How- 
ever, the actual impact of the tides on a single cluster-sized element may vary 
from the statistics and is generally much more complex. As an example, let's 
consider one of our stellar particle, orbiting in the merger. It will likely ex- 
perience a series of on-off compressive episodes, with different periods. On a 
speculative note, it is possible to imagine that external gas and dust would be 
driven by the compressive modes into the cluster core and then, condensed to 
form new stars. If the compressive events are well separated in time (and thus, 
in space along the orbit of the cluster), the chemical properties of this additional 
gas may differ from those of the original interstellar medium that formed the first 
generation of stars. In this fashion, it is possible to explain multiple populations, 
with diffe rent helium abuiidances, as observed in some resolved massive clusters 
(see e.g. iPiotto et al.ll2007l : iMilone et 311120091 ). Furthermore, it is possible that 
slow stellar ejecta would be kept bound by the additional gravitational effect of 
the compressive tides, and then be processed into one or more secondary gen- 
erations of stars. This would favor a self-enrichment of the cluster with respect 
to a tidally unperturbed system. Yet, a quantitative study of the impact of the 
compressive modes on the hydrodynamics of a star cluster is at a too early stage 
and this point remains speculative. 



6 Conclusions 



In this contribution, we showed that fully compressive tidal modes develop in 
the course of a prototypical merger like the Antennae galaxies. Accounting for 
~ 15% of the stellar mass, they last long enough (-^ 10 Myr) to prevent the 
newly formed cluster from early dissolution. A good correlation has been found 
between the position of the compressive regions and the observed substructures 
(star clusters and TDGs). When trying to extend these conclusions to a larger 
set of parameters, we showed that the fly-bys would not be strongly affected 
by the compressive modes. Indeed only mergers present a significant increase 
in the mass fraction entering such modes, matching the epochs of the orbital 
passages. Furthermore, the amplitude of these increases is anti-correlated with 
the pericenter distance: short distances induced a large rise in the mass fraction 
in compressive mode, while wide interactions are closer to the quiescent stage. 
Typical proto-clusters also show an on-off behavior of the compressive tides. 
This might give rise to multiple star formation epochs, with or without self- 
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enrichment of the cluster. However, the exact response of parsec-scale structures 
embedded in a tidal field varying in time and direction has still to be addressed 
and will be presented in a forthcoming paper. 
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